Gamow- Teller strength distributions at finite temperatures and electron capture in 

stellar environments 
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We propose a new method to calculate stellar weak-interaction rates. It is based on the Thermo- 
Field-Dynamics formalism and allows the calculation of the weak-interaction response of nuclei at 
finite temperatures. The thermal evolution of the GT+ distributions is presented for the sample 
nuclei 54 ' 56 Fe and 76,7S,80 Ge. For Ge we also calculate the strength distribution of first-forbidden 
transitions. We show that thermal effects shift the GT+ centroid to lower excitation energies 
and make possible negative- and low-energy transitions. In our model we demonstrate that the 
unblocking effect for GT+ transitions in neutron-rich nuclei is sensitive to increasing temperature. 
The results are used to calculate electron capture rates and are compared to those obtained from 
the shell model. 

PACS numbers: 26.50,+x, 23.40.-s, 21.60.Jz, 24.10.Pa 



I. INTRODUCTION 

The properties of nuclei at finite temperatures have 
attracted attention for a long time. Due to a consider- 
able amount of experimental data the main subject of the 
study was the thermal properties of the giant dipole reso- 
nance (see, e.g., and reference therein). In the astro- 
physical context the thermal properties of Gamow- Teller 
(GT) transitions are of special interest since they play a 
crucial role in weak-interaction mediated reactions (elec- 
tron capture, beta-decay, neutrino scattering etc.) [|[. 
For example, electron capture on iron group nuclei initi- 
ates the gravitational collapse of the core of a massive star 
triggering supernova explosions. Moreover, electron cap- 
ture rates largely determine the mass of the core and thus 
the fate of the shock wave formed by the supernova explo- 
sion. Since the strong phase space dependence makes the 
relevant stellar weak-interaction rates at the early stage 
of the collapse - when the electron chemical potential fi e 
is of the same order as the nuclear Q-values - very sensi- 
tive to the GT distributions these need to be calculated 
very accurately in this regime. With proceeding collapse 
and hence increasing density, fi e grows faster than the 
Q-values of the nuclei present in the matter composition, 
the capture rates become less sensitive to the details of 
the GT distribution and are mainly determined by the 
total GT strength and its centroid energy. However, for- 
bidden transitions can no longer be neglected when /i e 
reaches values of the order of 30 MeV at core densities 
p > 10 11 g/cm 3 (EH The situation is further complicated 
by the fact that weak-interaction processes in stellar en- 
vironments take place at temperatures of the order a few 
hundred keV to a few MeV and GT transitions occur not 
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only from the nuclear ground state, but also from excited 
states. 

From a microscopic point of view, there are two 
routes of handling GT strength distributions and weak- 
interaction rates at finite temperatures. One is 
a state-by-state evaluation of the rate by summing 
over Boltzmann-weighted, individually determined GT 
strengths for the various states. The second is based 
on an equilibrium statistical formulation of the nuclear 
many-body problem. In this approach, the thermal 
response of a nucleus to an external perturbing field 
is given by the canonical (or grand canonical) expec- 
tation value of the corresponding transition operator. 
When applied to charge-exchange processes this method 
yields the temperature-dependent GT and first forbid- 
den strength function that can be then used to calculate 
weak-interaction rates. 

For sd and p/-shell nuclei the first approach was origi- 
nally used by Fuller et al. @ who calculated stellar weak- 
interaction rates using the independent-particle shell- 
model, supplemented by experimental data, whenever 
available. To allow for GT transitions from nuclear ex- 
cited states these authors employed the Brink hypoth- 
esis. This assumes that the GT strength distribution 
on excited states is the same as for the ground state, 
only shifted by the excitation energy of the state. These 
rates were subsequently updated, taking into account 
the quenching of axial coupling constant Q. Modern 
high-performance computing capabilities combined with 
state-of-the-art diagonalization approaches make possi- 
ble shell-model calculations of the GT strength distribu- 
tion not only for the nuclear ground state, but also for 
the few lowest excited states. It was demonstrated 0, H| 
that even for the lowest excited states in the parent nu- 
cleus the Brink hypothesis is valid only for the bulk of 
the GT strength, but is not applicable for the individ- 
ual transitions to states at low-excitation energy in the 
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daughter nucleus. Taking this into account, the weak- 
interaction rates based on the shell-model diagonaliza- 
tion approach [3, [13] were derived from the individual 
GT distributions from the lowest excited states and from 
'back-resonant contributions' , i.e. from transitions de- 
termined from the inverse GT distributions connecting 
excited states in the daughter to the lowest states in the 
parent spectrum. However, the compilation of [3.[ic| ap- 
plied Brink's hypothesis when taking into account GT 
transitions from highly excited states. Weak-interaction 
rates have also been computed using the spectral dis- 
tribution theory 111! and the proton-neutron quasiparti- 
cle RPA model [12|. The first method is also based on 
the Brink hypothesis. The latter does not use this hy- 
pothesis, but some uncertainties arise due to the approx- 
imate treatment of the parent excited states as multi- 
quasiparticle states as well as the insufficient knowledge 
of the quantum numbers of the states involved. Recently 
stellar electron capture rates have been calculated within 
the framework of the finite temperature RPA using a set 
of Skyrme interactions 13]. This method has the ad- 
vantage of consistency, however, it misses relevant cor- 
relations which, as we will demonstrate in the present 
paper, are crucial to derive stellar electron capture rates 
for neutron-rich nuclei. 

The statistical way for the calculation of temperature- 
dependent GT and first-forbidden strength functions was 
first applied in [4] to study electron capture on neutron- 
rich nuclei. The most advances realization of this way is 
presently performed in the framework of the shell-model 
Monte-Carlo (SMMC) method 14]. It was found [H 
that with increasing temperature the GT centroids shift 
to lower excitation energies and the widths of the dis- 
tributions increase with the appearance of low-lying 
strength. Both effects arise from thermally excited states, 
i.e., the Brink hypothesis is not supported by SMMC cal- 
culations. In spite of its advantages, the SMMC method 
only yields the lowest moments of the GT strength distri- 
butions which introduce some inaccuracies into the rate 
calculations. Moreover, the SMMC method has restric- 
tions in its applicability to odd-odd and odd-A nuclei at 
low temperatures. 

Thus, the problem of an accurate description of the 
GT strength distribution at finite temperatures and re- 
liable estimates of stellar weak-interaction rates is not 
solved completely yet: The shell-model diagonalization 
approach allows for detailed spectroscopy, but partially 
employs the Brink hypothesis. The SMMC method is 
free from this disadvantage, but cannot provide a de- 
tailed strength distribution. Moreover, present computer 
capabilities allow the application of the shell-model di- 
agonalization method only to nuclei in the iron region 
(A = 45 — 65), whereas the SMMC approach can in 
principle be applied to all nuclei. However, such cal- 
culations are rather time-consuming and have therefore 
been limited to about 200 nuclei with mass numbers 
A = 65 — 120 [1, [Hj], while weak processes in more mas- 
sive and neutron-rich nuclei also play an important role 



in various astrophysical scenarios. Therefore, alternative 
methods to deal with Gamow- Teller strength distribu- 
tions and weak-interaction rates at finite temperatures 
are desirable. 

In this paper we study the temperature dependence of 
the Gamow- Teller strength applying the proton-neutron 
quasiparticle RPA [l7| extended to finite temperatur e by 
the Thermo-Field-Dynamics (TFD) formalism [18|, 
This technique has the advantage that it does not rely 
on Brink's hypothesis. The energies of the GT transi- 
tions and corresponding transition strengths are calcu- 
lated as functions of the nuclear temperature. In this pa- 
per, we apply this method for the calculations of weak- 
interaction rates on iron group nuclei and neutron-rich 
nuclei beyond pj-shell. However, this method is not only 
restricted to these nuclei but can be applied to heav- 
ier nuclei as well. In addition, it allows calculations of 
the strength distributions for forbidden transitions which 
contribute significantly to weak-interaction rates at high 
densities. Although, in the present paper, we restrict 
our study to the one-phonon approach, the method can 
be extended to higher phonon admixtures thus yielding 
more detailed strength distributions. 

The paper is organized as follows. In Sec. [Til some 
important features of the TFD formalism with applica- 
tion to the nuclear structure problem at finite tempera- 
tures are presented. The finite temperature (TQRPA) 
equations, which describe the strength distribution of 
charge-exchange transitions in hot nuclei, are given in 
this section as well. In Sec. IIIII the necessary formulae 
to calculate electron capture rates in a stellar environ- 
ment are introduced. Results for the GT strength dis- 
tributions and electron capture rates in 54 ' 56 Fe are pre- 
sented in Sec. lIVl Here, we also compare the results with 
those from the shell-model diagonalization approach. In 
Sec. [V] we study the temperature dependence of GT 
and first-forbidden strength distributions in the neutron- 
rich isotope 76 Ge. The corresponding electron capture 
cross sections and rates are calculated and compared 
with those obtained based on a hybrid SMMC+RPA 
model [20(. Conclusions are drawn in Sec. I VII 



II. FORMALISM 

As a method to study a thermal behavior of quan- 
tum many-body systems the TFD method has two at- 
tractive features: a) temperature effects arise explicitly 
as T-dependent vertices, providing a convenient start- 
ing point for various approximations; b) temperature and 
time are independent variables. The first feature allows 
for straightforward extensions of well-established zero- 
temperature approximations. It has been employed pre- 
viously in [111, [23, [^HH [2^] where selected nuclear struc- 
ture problems at finite temperatures were considered. 

The standard TFD formalism treats a many-particle 
system in thermal equilibrium with a heat bath and 
a particle reservoir in the grand canonical ensemble. 
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The thermal average of a given operator A is calcu- 
lated as the expectation value in a specially constructed, 
temperature-dependent state |0(T)), which is termed the 
thermal vacuum. This expectation value is equal to the 
usual grand canonical average of A. 

To construct the state |0(T)), a formal doubling of 
the system degrees of freedom is introduced. In TFD, 
a tilde conjugate operator A - acting in the independent 
Hilbert space - is associated with A, in accordance with 
properly formulated tilde conjugation rules [l8L fl9L f26j] . 
For a system governed by the Hamiltonian H the whole 
Hilbert space is now spanned by the direct product of 
the eigenstates of H (H\n) = E n \n)) and those of the 
tilde Hamiltonian H, both having the same eigenvalues 
(H\n) = E n \n)). In the doubled Hilbert space, the ther- 
mal vacuum is defined as the zero-energy eigenstate of the 
so-called thermal Hamiltonian Ti = H — H and it satis- 
fies the thermal state condition fijl USL \2&i 



A\0(T)) = ae n / 2T A*\0(T)), 



(1) 



where a = 1 for bosonic A and a — i for fermionic A. 

The important point is that in the doubled Hilbert 
space the time-translation operator is not the initial 
Hamiltonian H, but the thermal Hamiltonian Ti. This 
means that the excitations of the thermal system are ob- 
tained by the diagonalization of Ti. As it follows from the 
definition of Ti each of its eigenstates with positive en- 
ergy has a counterpart - the tilde-conjugate eigenstate - 
with negative energy, but the same absolute value. This 
is a way to treat excitation- and de-excitation processes 
at finite temperatures within TFD. 

Obviously, in most of practical cases one cannot di- 
agonalize Ti exactly. Usually, one resorts to certain ap- 
proximations such as the Hartree-Fock Bogoliubov mean 
field theory (HFB) and the Random-Phase Approxima- 
tion (RPA) (see e.g. [HI)- ^ n what follows the formal 
part of the TFD studies for charge-exchange excitations 
in hot nuclei is based in part on the results of [III, HI]. 

In our present study we use a phenomenological nu- 
clear Hamiltonian consisting of a static mean field, BCS 
pairing interactions and separable multipole and spin- 
multipolc particle- hole interactions, including isoscalar 
and isovector parts. This is usually referred to as the 
quasiparticle-phonon model (QPM) [29]. It was used 
to study the charge-exch ang e excitations in nuclei at 
zero temperature in [3(| \3l\ . In principle, the QPM 
formalism enables one to go beyond the QRPA and 
take into account the coupling of quasiparticles and 
phonons. At finite temperatures this coupling was con- 
sidered in [24|, [13] • However, in the present study we 
restrict ourselves to the thermal QRPA. 

The main line of the present discussion is very similar 
to the QPM at T = ||. We begin with the thermal 
Hamiltonian, which reads as 



Hqpm — Hqpm — Hqpm — 7~Ls 



n 



(2) 



The first step in the approximate diagonalization of 
7Yqpm is the treatment of the pairing correlations. This 



is done by two successive unitary transformations. The 
first is the usual Bogoliubov u, v transformation from the 
original particle operators a] m , a lm to the quasiparticle 



ones a 



ol-. The same transformation is applied to 



jm 

tilde operators aj m , a,j m , thus producing the tilde quasi- 
particle operators a\ m , 3 jm . 

The second transformation is the so-called thermal Bo- 
goliubov transformation [l8l . [l9| . It mixes the quasipar- 
ticle and tilde quasiparticle operators, thus producing 
thermal quasiparticle operators and their tilde partners: 

4m> #/m> Pjmi Pjm' We use Ojima's [1| complex form 
01 the thermal Bogoliubov transformation 



@}m = x J a l 



x 3 01 jm 



(*3 



(3) 



The reasons for this are given in [271 ]. 

The coefficients Uj, Vj, Xj, yj are found by diago- 
nalizing Ti sp + 7i pa ir and demanding that the vacuum 
of thermal quasiparticle obeys the thermal state condi- 
tion ([I]). This is equivalent to the minimization of the 
thermodynamic potential for Bogoliubov quasiparticles. 
As a result one obtains the following equations for u 
and Xj, yj\ 



3 1 



1 

71 



Ei — X T 



1/2 



(1 



Vi 



1 + cxp 



-1/2 



2\l/2 



(4) 
(5) 



where Sj = y/ (Ej — A T ) 2 + and r is the isospin quan- 
tum number r = n,p. 

The pairing gap A r and the chemical potential A T are 
the solutions to the finite-temperature BCS equations 

A.(T) = ^£ T (2j+l)(l-2y>^, 



E T (2 J + 1)(^+ U ^), 



(6) 



where N T is the number of neutrons or protons in a nu- 
cleus and ~Y^ T implies a summation over neutron or pro- 
ton single-particle states only. From the numerical solu- 
tion of Eqs. ([6]) it is found that the (pseudo)critical tem- 
perature is T cr w iA T (0) (see e.g. [32|,[33j]) in accordance 
with the BCS theory. 

With the coefficients Uj, Vj, Xj, yj, defined by 
Eqs. dU E]), the one-body part of H sp + 7i pa ir reads 

H sp + W pair ~ £ ]T T £j (0} m f3 jm - ]3]J jm ) (7) 

r jm 

and corresponds to a system of non-interacting thermal 
quasiparticles. The vacuum for thermal quasiparticles 
(hereafter denoted by |0(T);qp)) is the thermal vacuum 
in the BCS approximation. The states /3j m |0(T); qp) 
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have positive excitation energies whereas the correspond- 
ing tilde-states /3j m |0(T); qp) have negative energies. 

The coefficients y| defined through (O determine the 
average number of thermally excited Bogoliubov quasi- 
particles in the BCS thermal vacuum 

(0(T);qp|aj m a im |0(T);qp) = y) (8) 

and, thus, coincide with the thermal occupation factors 
of the Fermi-Dirac statistics. Since the thermal vacuum 
|0(T); qp) contains a certain number of Bogoliubov quasi- 
particles, excited states can be built on |0(T);qp) by ei- 
ther adding or removing a Bogoliubov quasiparticle. Be- 
cause of 

a]ml°( T );qp> = ^/3jml°( T );qp); 

a^|0(T);qp) = -%^|0(T);qp) 

(a m =(-l) j - m ^J (9) 

the first process corresponds to the creation of a non-tilde 
thermal quasiparticle with positive energy, whereas the 
second process creates a tilde quasiparticle with negative 
energy. 

In the next step of the approximate diagonalization of 
T^qpm, long-range correlations due to the particle- hole 
interaction are taken into account within the proton- 
neutron TQRPA. The part of H p h in ^ responsible for 
charge-exchange excitations reads 

H% = -2 4 A) (<,M A/I - MlM X)i ) 

-2j24 LX) (Sl Xti S LXti - Sl Xfl S LXfi ) , (10) 

LA/i 

where M Xfi and Sj^ x are single-particle multipole and 
spin-multipole operators: 

j p m p 

S L\» = E (ip m ?l iArA [ F i°']M t(_) l-?™ TO ™) im p Snrn„ ! 
j p m p 

[Y L <r]* = ^(LMlmlA^yiM^,^)^. (11) 

M, m 

The parameters k\ and denote the strength pa- 

rameters of the isovector multipole and spin-multipole 



forces, respectively. The states of natural parity are gen- 
erated by the multipole and spin-multipole L = A interac- 
tions, while the spin-multipole interactions with L = Ail 
are responsible for the states of unnatural parity. 

Within the TFD formalism the TQRPA equations are 
derived in the following way. First, H p \ is written in 
terms of the thermal quasiparticle operators. Then, 
the sum of ([7]) and HQ is diagonalized with respect to 
charge-exchange thermal phonons. 

The thermal phonon creation operator is defined 
as a linear superposition of the proton-neutron thermal 
two-quasiparticle operators 

+ Ifc&nft + > (12) 

where [ denotes the coupling of single-particle angu- 
lar momenta jmjp to total angular momentum A. Now 
the thermal equilibrium state is treated as the vacuum 
|0(T); ph) for the thermal phonon annihilation operators. 

The thermal phonon operators are as- 
sumed to commute as bosonic operators, i.e., 

[Qx^iiQ\'u'i'] = £u"W<fe'- Tnis assumption imposes 
the following constraint on the phonon amplitudes: 

E (v&.i& + + <t, 'w 

jpjn 

Furthermore, the phonon amplitudes obey the closure 
relation. 

Demanding that the vacuum of thermal phonons obeys 
the thermal state condition ([T]) and applying the varia- 
tional principle to the average value of thermal Hamilto- 
nian with respect to one-phonon states Qw|0(r); ph) 

or Q^_|0(T);ph)[Il] under the constraints |(T3J) one 

gets a system of linear equations for the amplitudes 

^i'in' ^j P j n > Vjpjn' ^ipj n ' etc - The system has a nontriv- 
ial solution if the energy u>xi of the thermal one-phonon 
state obeys the following secular equation: 
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A' 
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(+-) y( 



(-+) 



X 
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(+-) 



X 



(-) 



-<H& ( 6 ) ^bb 



( + -) 



(-+) v( 



A 66 _ 



.(6) 



o, 



(14) 



r 



where a = A and b = XX for excitations of natural parity, 
while a = (A — 1)A and 6 = (A + 1)A for unnatural parity 



excitations. The functions X 
d = a, b) in (fl4|) are defined as 



(±) 

cd ' 



X 



(±=F) ^_ 



(c = a, b and 



9 f F i+) (V {±) ) 2 

Y (± )/ N _ i V"^ f ( C ) ,(d) I £ j P jA U j P j n ) 

x cd W - -jZ^ • / ., .- ./ 1 77+) 



) 2 -CJ 2 



|2 — ,,i2 



v(±t) / ,\ _ 2w V f (c) f (d) / 
j v 3n y 



JpJn 

2 



-51 -51 

JpJn JpJn 

[ef] ) 2 - c^ 2 



L. (15) 



Here f>*\- and f,)^ denote the reduced single- 

" JpJn, ** Jp Jti u 



(AL) 



particle matrix elements of the multipole and spin- 
multipole operators (fTTj) ; Ujj n = u j P v j n ± v j P u j n an d 

(±) - - ■ -<- - ■ • ^ - V2A + 1. 



(±) _i_ 

JpJn JP J" JP J 



£ jp ± £ 3~ ; ^ 



Let us consider the secular equation in detail. The 
poles ej; which do not exist in the QRPA equations 

J P Jn 

at zero temperature arise from the crossed terms /3^/3^ 
in the thermal phonon operator definition (|12|l . Due to 
these poles, new states appear in a low-energy part of 
the thermal excitation spectrum. In contrast to the zero 
temperature case, the negative solutions of the secular 
equation now have a physical meaning. They correspond 
to the tilde thermal one-phonon states and arise from 
[fl terms in the thermal phonon operator. As it was 
noted above, creation of a tilde thermal quasiparticle cor- 
responds to the annihilation of a thermally excited Bo- 
goliubov quasiparticle. Consequently, excitations of low- 
and negative-energy thermal phonons correspond to tran- 
sitions from thermally excited nuclear states. Further- 
more, when the pairing correlations vanish (i.e., T > T cr ), 
some poles no longer contribute to the secular equation 
since the corresponding numerators vanish. This is true 
for particle-particle and hole-hole Sjj n poles as well as 



for particle-hole e\ } poles 

JpJn 



The expressions for the thermal charge- exchange 
phonon amplitudes can be found in (28j. The ampli- 
tudes depend on both the quasiparticle and the phonon 
thermal occupation numbers. Some remarks are in order. 
TQRPA equations for GT excitations at finite tempera- 
ture were also derived in [34j . They were obtained using 
the equation of motion method by replacing vacuum ex- 
pectation values by thermal averages, i.e., without apply- 
ing the TFD formalism and doubling the Hilbert space. 
Therefore, in contrast with the present study negative so- 
lutions of the TQRPA equations were neglected in [lij . 

After diagonalization within the TQRPA the thermal 
Hamiltonian 7Yqpm becomes 



Hqpm 



). (16) 



The vacuum |0(T); ph) of thermal phonons is the thermal 
vacuum in the thermal quasiparticle RPA. Since we use 
the thermal BCS approximation, which violates the par- 
ticle number, the charge- exchange thermal one-phonon 
states are superpositions of states, which belong to the 
daughter nuclei (N-l, Z+l) and (N+l, Z-l). They de- 
couple at temperatures T > T cr , when the pairing corre- 
lations vanish. Then, if the state Q\ -|0(T);ph) belongs 
to the (N ± 1, Z T 1) nucleus, the state Ql_|0(T); ph) is 
in the (N =p 1, Z ± 1) nucleus. 



III. ELECTRON CAPTURE RATES 

Considering electron capture in stellar environments 
we make the following assumptions: 1. the temperature 
in the stellar interior is so high that atoms are fully ion- 
ized, and the surrounding electron gas is described by a 
Fermi-Dirac distribution, with temperature T and chem- 
ical potential /i e . Neutrinos escape freely from the inte- 
rior of the star. Hence no Pauli blocking for neutrinos 
is considered in the final state. 2. the parent nucleus 
is in a thermal equilibrium state treated as the thermal 
(phonon) vacuum. 3. electron capture leads to charge- 
exchange transitions from the thermal vacuum to ther- 
mal one-phonon states. 

Under these circumstances the electron capture rate is 
the sum of the transition rates from the thermal vacuum 
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to the i-th thermal one-phonon state of the multipolar- 
ity J 



A c 



In 2 



6150 sec 



EE<^ 0C = EE^ 



(17) 



j % 



J i 



Here $ is the squared reduced matrix element of the 
transition operator between the thermal phonon vacuum 
and a thermal one-phonon state (see below); F° c is a 
phase space factor which depends on the transition en- 
ergy -Ej]~ and can be found elsewhere @. 

Denoting the proton-to-neutron (p — > n) transition op- 
erator with multipolarity J as Dj one obtains the fol- 



lowing expression for the transition strength $ 



(+) 

Ji 



3> 



( + ) 



,,, - | (0(r);ph||Q. 7A « J DS +) H0(r);ph) 



3p3n 



where dy\j p j n ) = (in II -Dj II jp) is a reduced single- 
particle matrix element of the transition operator, and 
the function £l(j p j n ; Ji) is given by 

n(jpjn] Ji) = v jp u jn {x Jp x ]n ip J jX + !l, !lj o'/;, ) 

+ u h u Jn (y jp x jn rj£ jn + x jpVj j£ jn ). (19) 

The transition strength to the tilde one-phonon state can 
be easily obtained from fT5|) and (|19p by changing non- 
tilde phonon amplitudes by their tilde counterparts and 
vise versa. (The expressions for the transition strengths 

$ j - \ corresponding to inverse n — > p transitions are given 
in [28] . It has been proved in 28] that the approach used 
here fulfills the Ikeda sum rule for Fermi and Gamow- 



p(+) 



Teller transitions.) 

The transition energy (parent excitation energy) E K U 
can be obtained from the energy shift between the proton 
subsystem of the parent nucleus and the neutron subsys- 
tem of the daughter nucleus including the proton-neutron 
mass difference. Thus we have 



E { jf = uj, h + (A/i 



lip 



Ar 



'tip 



), 



(20) 



where A/i„ p = /i p — \i n is the difference between neutron 
and proton chemical potentials and Am np = m n — m p 
is the neutron-proton mass splitting. Note, that at finite 
temperature the energies as well as ojji can be both 
positive and negative. Thus, to the capture process ther- 
mal one-phonon states with both positive and negative 
values of E^ contribute to the rate. 

In what follows in (jTTJ) we take into account the con- 
tributions from allowed (Gamow- Teller and Fermi) tran- 
sitions and first- forbidden transitions. The operators of 



allowed Fermi and Gamow- Teller transitions are taken in 
the standard form 



D$=g A <rt<-+\ 



(21) 



where is the isospin raising operator. For the oper- 
ators of the first-forbidden n — ► p transitions the non- 
relativistic form is used 



rr ■ p aZ 

-i 1 

2R 



m 



(+) 



D 



D 



(+) = 
1- 

(+) 



p aZ . . 

9v 7T5\9av x r - ig v r) 

m 2K 



■9A r 
\/5 



(22) 



In Eqs. (|2"Tj) and (|22[) r, p, and er refer to the coordinate, 
momentum and spin operators of a nucleon, gy — 1 and 
g A = —1.25 denote the vector and axial coupling con- 
stants, a is the fine structure constant, R are the 
charge and the radius of the nucleus; m is the nucleon 
mass and p e and q v denote the momenta of the incoming 
electron and outgoing neutrino, respectively. 



IV. IRON ISOTOPES 

In this section we discuss the numerical results for the 
p/-shell nuclei 54 ' 56 Fe. Experimental data are available 
for these nuclei to test our calculations at zero tempera- 
ture. Moreover, these iron isotopes are among the most 
essential nuclei in their importance for the electron cap- 
ture process for the early presupernova collapse @, [H| . 

The proton and neutron mean fields are described 
by spherically symmetric Woods-Saxon potentials with 
parameters from (36|. We only readjust the poten- 
tial depths to fit the proton and neutron binding en- 
ergies of the parent nucleus to their experimental val- 
ues. The single-particle basis includes all discrete 
bound states as well as selected quasi-bound states 
with large j in the continuum. The proton (neutron) 
pairing strength parameters G p t n ) are fixed to repro- 
duce the odd-even mass difference through a four term 
formula [3^ ] involving the experimental binding ener- 
gies (38|. At T = the obtained proton and neu- 
tron BCS energy gaps are: A p ( n ) = 1.52 (0.0) MeV for 
54 Fe,[l| A p{n} = 1.57 (1.36) MeV for 56 Fe. The isovec- 

tor strength parameters k^ 01 ^ and k^ 21 ^ are adjusted 
to reproduce the experimental centroid energies of the 
GT_ and GT + resonances in the nuclei under consider- 
ation [39l. |40|. [4l| . The corresponding values of k^ 01 -* and 

k± are in agreement with the rough estimates in [42| . 
The spin-quadrupole interaction weakly affects the GT 
strength distributions. 

The total GT strengths calculated with the bare GT± 
operators at^ are S+ = 6.6, S- = 12.7 in 54 Fe and 
5+ = 5.1, S- = 17.0 in 56 Fe respectively. These S T 
values obey the Ikeda sum rule (a small deviation is 
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FIG. 1: (Color online) Left panels: Comparison of GT+ ex- 
perimental data pp| . l4lj with the calculated QRPA strength 
distributions for 54 ' 56 Fe. The QRPA peaks are scaled by 0.5 
for convenience. Right panels: Comparison of the GT+ run- 
ning sums corresponding to the experimental, QRPA, and 
LSSM [{| strength distributions. 



due to the incompleteness of our single particle basis) 
but noticeably overestimate experimental data (see e.g. 
[HE S3)- This is common for any RPA (or QRPA) 
calculation of the GT strength and is remedied by an ef- 
fective value for the axial weak coupling constant. We 
use g* A = Q.7AgA as in shell-model calculations @. 

In FigJTJ the experimental and theoretical (quenched) 
distributions of GT + strengths are presented. Here we 
also compare the GT + running sums corresponding to 
the experimental, QRPA, and large scale shell-model 
(LSSM) strength distributions. One can see that 
the QRPA calculations reproduce the resonance positions 
but not the fragmentation of the strength. It is a well 
known fact that RPA calculations cannot describe the full 
resonance width (at least in spherical nuclei) and produce 
only a part of it, the so-called Landau width. The latter 
is quite small for the GT resonance. As a result, the near 
threshold part of the GT + strength, which corresponds 
to transitions to low- lying 1 + states in the daughter nu- 
clei 54,56 Mn, is not reproduced in our calculations. In 
this respect the shell-model calculations are clearly at an 
advantage. 

We now turn to the temperature evolution of the GT + 
strength distributions. The strength distributions for 
54,56p e a -j- severa i temperatures are shown in Fig. O All 



figures are plotted as a function of the energy transfer E 
to the parent nucleus. 

With increasing temperature, in our model, two effects 
occur that influence the GT + strength distribution: 

(i) At low temperatures, due to pairing, GT + transi- 
tions involve the breaking of a proton Cooper pair associ- 
ated with some energy cost. This extra energy is removed 
at T > T„ (T cr 0.8 MeV) and the peak in the GT+ 
distribution moves to smaller energies. Some extra en- 
ergy has to be paid at low temperatures to add one more 
nucleon to the neutron subsystem of 56 Fe because of a 
non-zero neutron energy gap. Obviously, this energy is 
also removed at T > T CT . 

(ii) GT_|_ transitions, which are Pauli blocked at low 
temperatures due to closed neutron subshells (e.g., I/7/2 
orbital), become thermally unblocked with increasing 
temperature. Similarly, protons which are thermally ex- 
cited to higher orbitals can undergo GT + transitions. In 
TFD such transitions are taken into account by fij 0L-, 

0j^f3j n , and Pj^0j^ components of the thermal phonon. 
Because of thermally unblocked transitions, some GT + 
strength appears well below the zero-temperature thresh- 
old, including negative energies. 

Due to the vanishing of the pairing correlations and 
appearance of negative- and low-energy transitions, the 
ccntroids of the GT + strength distributions in 54 > 56 Fe are 
shifted to lower excitation energies at high temperatures. 
Our calculations indicate that a temperature increase to 
0.8 MeV results in the GT + centroid shifts of the or- 
der of 1.5 MeV for 54 Fe and 2.5 MeV for 56 Fe. Thus 
the present approach violates Brink's hypothesis. Simi- 
lar results have been obtained in SMMC calculations of 
the GT centroids at finite temperatures. We also ob- 
serve (see Fig. [2]) a gradual decrease of the total GT + 
strength when the temperature increases from zero to 
0.8 MeV. Nevertheless, as was pointed out above, the 
present approach preserves the Ikeda sum rule at finite 
temperatures. 

The calculated GT + strength distributions have been 
used to obtain the stellar electron capture rates for 
54 ' 56 Fe. The rates have been calculated for densities 
between log(pV e ) = 7 and \og(pY e ) = 10 as a function 
of temperature T 9 (T 9 = 10 9 K and 1 MeV« 11.6 T 9 ). 
The comparison between the TQRPA rates and the large- 
scale shell-model results [Io| is presented in Fig. [3] 

As it must be, the electron capture rates increase with 
temperature and density. Due to the larger value of the 
zero-temperature threshold Q for 56 Fe, both approaches 
yield a higher rate for 54 Fe than for 56 Fe at a given tem- 
perature and density. Both approaches give very similar 
values for the strength and the location of the GT + res- 



onance m 



54.50 



Fe at T — 0. Therefore, the excellent 
agreement between the TQRPA and shell model rates 
at log(pV e ) = 10 and low temperatures (/i e sa 11 MeV) 
is not surprising, since the rates are dominated by the 
resonance contribution. 

The more interesting point is that at high temperatures 
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FIG. 2: (Color online) Temperature evolution of GT+ strength distributions for 54 Fe (upper panels) and for 56 Fe (lower panels) 
versus parent excitation energy. The full (dashed) lines refer to transitions to non-tilde (tilde) thermal one-phonon states. St 
is the total GT+ strength. The arrows indicate the zero-temperature threshold Q = Mf — Mi, where Mij are the masses of 
the parent and daughter nuclei. Q( 54 Fe) = 1.21 MeV and Q( 56 Fe) = 4.20 MeV. 



the TQRPA rates always surpass the shell-model ones. 
To understand this point, it needs to be clarified which 
part of the TQRPA GT + strength dominates the electron 
capture at a given temperature and density. To this end 
we calculate the relative contributions A° c /A cc of differ- 
ent thermal one-phonon states to the capture rates for se- 
lected values of temperature and density, (Tg,\og(pY e )). 
The results are depicted in Fig. |U 

At low temperatures and densities (Fig. |UJa),(b)), i.e., 
when p e is small and high-energy electrons from the tail 
of the Fermi-Dirac distribution are not sufficiently avail- 
able to allow for efficient capture on the GT + resonance, 
the TQRPA capture rates are dominated either by the 
negative-energy ( 54 Fe) or by the low-energy ( 56 Fe) part 
of the GT + strength which originates from thermally 
unblocked p — > n transitions. The TQRPA rates are 
larger than those of the shell-model at low (T, p) due 
to differences in the strength and the energy of such 
transitions. We note that in the shell model evalua- 
tion negative-energy transitions were mainly included by 
back-resonances; i.e. by inverting the Fermi and GT_ 
strength distribution of 54 Mn and 56 Mn, respectively. 
Different to the TQRPA approach the shell model GT_ 
distributions of these nuclei are highly fragmented due to 



correlations and have centroids at rather high excitation 
energies in 54 Fe and 56 Fe which, at low temperatures, are 
strongly suppressed by the Boltzmann factor. In partic- 
ular, the differences in energy positions of the transitions 
are important since at low (T, p) the rates can change 
drastically by a small change in a transition energy. To 
see whether the TQRPA reliably predicts the energy and 
the strength of negative- and low-energy transitions one 
needs to go beyond the TQRPA. 

At \og{pY e ) = 9 and T 9 < 5 (p e » 5.1 MeV) the 
near threshold part of the GT + strength dominates the 
capture rates. Since this part is not reproduced within 
the TQRPA, the rates appear to be smaller than the 
LSSM ones. To test this hypothesis, the capture rates 
at log(pY~ e ) = 9 have been calculated - guided by the 
shell-model GT + distributions @ - assuming that the 
near threshold GT+ strengths for 54 Fe and 56 Fe are 0.1 
and 0.2, respectively. We therefore assign the value 0.1 
(0.2) to the GT+ strength in 54 Fe ( 56 Fe) at the zero- 
temperature threshold. (A similar method was used 
in p, @ to include the contribution of low-lying tran- 
sitions.) This yields a much better agreement between 
the TQRPA and shell- model rates (see Fig. [3]). Thus, 
to improve the reliability of the TQRPA for the cal- 
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FIG. 3: (Color online) Electron capture rates for 54 Fe (upper panels) and 56 Fe (lower panels) calculated using the TQRPA 
approach as a function of temperature (Tg measures the temperature in 10 9 K) and for selected values of density pY e (in 
gem -3 ). For comparison, the LSSM rates are also shown. The dashed lines at \og(pY e ) = 9 correspond to the TQRPA 
rates calculated with the assigned near threshold strength (see text). 



culation of stellar electron capture rates at intermedi- 
ate densities and low temperatures, the fragmentation of 
the GT + resonance should be considered to reproduce 
the near threshold GT + strength. At higher tempera- 
tures the near threshold strength becomes less impor- 
tant. In Fig. d£c),(d) the relative contributions A? c /A ec 
at (Tg, log(pY" e )) = (5,9) with and without the assigned 
near threshold strength are depicted. As can be seen, 
the contribution from the near threshold strength is not 
dominant. 

When the temperature approaches Xg ps 10, the 
rates are dominated by the strong transitions involv- 
ing the GT+ resonance at low (Fig. [He)) as well as at 
high (Fig.EKf)) densities. (Note that at \og(pY e ) = 7 and 
hight temperatures the contribution of negative-energy 
transitions is non- negligible.) As was discussed above, 
the TQRPA predicts that with increasing temperature 
the GT + resonance shifts to lower excitation energies. 
This explains why the TQRPA rates always surpass the 
LSSM ones at high temperatures. 

Thus, only at high densities and low temperatures the 
TQRPA and LSSM electron capture rates for 54 > 56 Fe are 
in a good agreement. As was mentioned above the dis- 
agreement at moderate densities and low temperatures 
can be removed by considering the fragmentation of the 
GT + strength. For a separable residual interaction used 
here this can be done following the method developed 



within the quasiparticle-phonon nuclear model, i.e., by 
taking into account the phonon coupling. The interesting 
question is how the phonon coupling affects the negative- 
and low-energy part of the GT + strength, which domi- 
nates the capture rate at low temperatures and densities. 
This is an open question and requires further investiga- 
tions. 



V. NEUTRON-RICH GERMANIUM ISOTOPES 

During gravitational collapse the nuclear composition 
moves towards higher mass number and more neutron- 
rich nuclei. Eventually nuclei will have all neutron pf- 
shell orbits filled, with valence neutrons in the sdg- 
shell (N > 40) and valence protons within the p/-shell 
(Z < 40). The Pauli principle blocks GT + transitions in 
such neutron-rich nuclei if the independent particle model 
is used. It has been demonstrated in [4( that at high 
enough temperatures, T ~ 1.5 MeV, GT + transitions 
become unblocked by thermal excitations which either 
move protons into the I59/2 orbital or remove neutrons 
from the p/-shell. An alternative unblocking mechanism, 
configuration mixing induced by the residual interaction, 
was considered in [20]. Based on this approach it was 
found that electron capture on nuclei with N > 40 is 
also dominated by GT + transitions even at rather low 
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FIG. 4: (Color online) Relative contribution \^ c /\ cc of z- 
th thermal one-phonon state to the electron capture rate 
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tions of the assigned near threshold strengths (see text). 



stellar temperatures, T ~ 0.5 MeV. Contrary to Q, it 
was argued that unblocking effects due to mixing are not 
too sensitive to increasing temperature. 

Consistent calculations of the electron capture rates 
for neutron-rich nuclei are not yet feasible in the shell 
model due to the large model space. In [20j the cap- 
ture rates have been calculated adopting a hybrid model: 
The partial occupation numbers calculated within the 
SMMC approach at finite temperature were used in cal- 
culations based on the RPA. In this subsection, using the 



germanium isotopes 



76,78,80 



Ge as examples, we apply the 



TQRPA formalism to calculate electron capture rates on 
neutron-rich nuclei. Particular attention is paid to the 
temperature dependence of the unblocking effect. 

For our TQRPA calculations the parameters of the 
model Hamiltonian for 76 > 78 > 80 Ge are chosen in the same 
manner as for 54,56 Fe. The sequence of single-particle 
levels obtained is close to that used in [H for 82 Ge. 
For pairing gaps we obtain: A p ( n ) = 1.50 (1.57) MeV 
for 76 Ge, Ap (n ) = 1.59 (1.42) MeV for 78 Ge, and 
A p(n) = 1.39 (1.35) MeV for 80 Ge. Our calculations take 
into account both the allowed (GT and Fermi) and first- 
forbidden transitions with J < 2. To generate one- 
phonon states of natural and unnatural parity we use 
the isovector multipole and spin-multipole strength pa- 
rameters K 



K-i 



(A = 0, 1, 2) according to [42j,|4S 



As a representative example, the strength distribution 
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FIG. 5: (Color online) Strength distribution (folded) of al- 
lowed (0 + and 1 + ) p — » n transitions in 76 Ge at various tem- 
peratures T. E denotes the transition energy. The contri- 
bution of + transitions is shown by the dashed line. St is 
the total strength. The arrows indicate the zero temperature 
threshold Q( 76 Ge) = M } - Mi = 7.52 MeV. The letters label 
the transitions: A = l/f /2 -» l/ 5 n /2 , B = \g\ j% -> lg$ /: 



^9/2 



7/2- 



of allowed p — » n transitions from 76 Ge for different val- 
ues of temperature is displayed in Fig. [5] The distri- 
butions have been folded with a Breit-Wigner function 
of 1 MeV width. As it follows from our study as well as 
from [3, [2(|, two single-particle transitions mainly con- 
tribute to the total GT + strength in neutron-rich ger- 
manium isotopes. These are the l<?g/ 2 ~~ > I37/2 particle- 



particle and I/7/2 ~* I/5/2 hole-hole transitions. In an 
independent particle model both transitions are blocked 
at zero temperature. However, in the present model they 
become unblocked due to pairing correlations and ther- 
mal excitations. Referring to Fig. [5] it is observed that, 
with increasing temperature, the peaks in the GT + dis- 
tribution shift to lower excitation energies and the total 
strength decreases in the vicinity of the critical temper- 
ature (T cr - 0.8 MeV) [Hoi. The shift is of the order of 
8 MeV and, hence, cannot be explained solely by remov- 
ing the extra energy needed to break a proton pair. 

To explain both effects we neglect the residual particle- 
hole interaction and consider the pairing interaction only. 
(As it follows from our study, the position of the GT + 
peaks in 76 > 78 ' 80 Ge is little affected by the inclusion of 
QRPA correlations and thermal one-phonon states can 
be considered as thermal two-quasiparticle states.) At fi- 
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nite temperature the GT + operator can excite configura- 
tions of four different types, namely, [flj^Pj^lln 
\fi\ 0j^^ , and [/jL/jjji . The respective transition ener- 
gies and transition strengths are 

Ei(j P -> jn) =£j p + e jn + Q*, 

*i& -> Jn) = {f} P u¥ v %AAA 

E2(j P ^jn)=-{Sj p +£ jn ) + Q*, 



E z{jp jn) 



'Op 



Q*, 



MJp - Jn) = (fU^lvlxlyl; 

Ei{jp -> Jn) = - (£j p - EjJ + Q*, 

Mj P -» jn) - (/LJ 2 «yl^„' 



(23) 



where Q* 
lowing j p 



Afi np + Am„p (see Eq. (j2"01) ), In the fol- 

l<7?/ 2 or the 



I/7/2 



1/: 



5/2 



j n refer either to the lc/g 
transition 



9/2 



At T < T cr the excitation of the /3j configuration 
dominates the strength distribution because of the fac- 
tor x 2 p x 2 n ~ 1 in $i(jp — * j n )- Therefore, at relatively 
low temperatures, when configuration mixing induced 
by pairing correlations in the ground state is the main 
unblocking mechanism, the position of the GT + peaks 
is given by #i(lff£ /2 -> lgy 2 ) and £i(l/£ /2 -> l/£) 2 ) 
(Fig. [S^a),(b)). With increasing temperature states hav- 
ing internal configurations other than those of the nuclear 
ground state gain statistical weight and, in particular, 
the pairing correlations in these excited states decrease. 
When the pairing correlations disappear and the factors 
v 2 p u 2 n in $i(jp — > j n ) become zero, the peaks considered 
completely vanish (Fig. (He)). The value of ^(jp — * jn) 
becomes zero as well. At T > T cr the poles £i g p /2 +£i 9 ™ /2 
and £i/p /2 + £ i/„" /2 n0 longer contribute to the secular 
equation HH). 

At T > T CI GT + transitions from (to) thermally oc- 
cupied (unblocked) orbitals dominate the strength distri- 
bution. Such transitions correspond to the excitation of 
the Id] P 0\„ n ]}, and \B\ tP d] fn ]}, configurations and 



9/2 



y 7 /2 J 



J 7/2 



their energies are E^(lg^^ 2 



1.9 7 " /2 ) and E 3 (lfP /2 - 
l/^, 2 ), respectively. Neglecting Am np , these energies are 
the energy difference between the final and initial single- 
particle states, i.e., E4(lg^ 2 — ► l<?"/ 2 ) ~ E\ g n._ — E 

and£ 3 (l/ 7 P /2 ->l/ ; 



9/2 



'5/2/ 



E 



Ei fp . Because of the 

J T/2 



thermally unblocked transitions the GT + peaks appear 
near the zero temperature threshold (Fig. [5Jd)). 

Thus, in contrast to [2(J, we find that the unblock- 
ing effect for GT + transitions in neutron-rich nuclei is 
sensitive to increasing temperature. No shift to lower ex- 
citation energies for the GT + peaks or decrease of the 
total GT + strength in the vicinity of the critical temper- 
ature were observed in 20]. To understand these differ- 
ences we compare again the approximations underlying 
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FIG. 6: (Color online) Upper panels: Occupation numbers for 
the I/7/2, I99/2 proton orbitals in 76 Ge as a function of tem- 
perature. Middle panels: Occupation numbers for the I/5/2, 
l<?7/2 neutron orbitals. Lower panels: The unblocking proba- 
bilities for the 1/7/2 — > 1/5/2 an d l<?g/ 2 ~~ * lfl™/2 transitions. 



the present model and the hybrid approach of [20] . The 
TQRPA has the virtue of consistency. It describes cor- 
relations by configuration mixing derived from a pairing 
interaction up to the 2p2h level. In the hybrid model 
occupation numbers at finite temperature are calculated 
within the SMMC approach, accounting for all many- 
body npnh correlations induced by a pairing+quadrupole 
residual interaction. These occupation numbers have 
then been used to define a thermal ground state which 
is the basis of an RPA approach to calculate the cap- 
ture cross sections, considering only lplh excitations on 
the top of this ground state. Therefore the hybrid model 
does not include explicitly 2p2h pairing correlations when 
calculating strength distributions. 

Repeating our observation from above, the TQRPA 
has two distinct transitions to overcome Pauli blocking. 
Using, for the sake of simplicity the language of the Inde- 
pendent Particle Model, in the TQRPA GT transitions 
can occur from a configuration-mixed states with OpOh 
and 2p2h components. These transitions lead to excited 
states with centroids which are shifted by the excitation 
energy of two particles which are raised across the pf—sdg 
shell gap, which corresponds to about 8 MeV for 76 Ge. 
The two peaks observed in the TQRPA GT strength dis- 
tribution in Fig. [5ja) ,(b) correspond to these two tran- 
sitions. As the 2p2h component has two neutron holes, 
GT transitions into these holes are not Pauli blocked. 
Hence these transitions are relatively strong within the 
TQRPA model at low temperatures. On the other hand, 
GT transitions between pure OpOh components are Pauli 
blocked. Transitions to final states corresponding to the 
lower centroid are only possible due to the small mixing 
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of 2p2h configurations into the final states. Hence the GT 
strength corresponding to the lower peak is rather weak 
at low temperatures (Fig. EJb)). The relative weight of 
the transition strength between these two peaks changes 
with increasing temperature due to the growing of the 
thermal excitations and the decreasing correlations in- 
duced by pairing. The later effect dominates at modest 
temperatures. As a consequence, the strength in the up- 
per peak decreases, while the lower peak grows and, at 
temperatures beyond the critical temperature T cr , domi- 
nate the GT strength distribution in the TQRPA model 
(Fig.Hd)). 

It is found that in the SMMC approach many-body 
correlations lead to much stronger excitation of particles 
(mainly neutrons) across the pf-sdg shell gap than found 
in the TQRPA model. This is demonstrated in Fig.^a)- 
(d), which compares the SMMC and thermal BCS occu- 
pation numbers for various orbitals as function of tem- 
perature. While BCS predicts only 0.3 neutrons to be 
excited out of the p/-shell at T = for 76 Ge, this number 
is about 1.2 for the SMMC, which is actually also smaller 
than the number of neutron excitations across the shell 
gap, derived recently experimentally (2.48 ± 0.30) [44| . 
Obviously the differences in occupation numbers leads 
to a larger Pauli unblocking within the SMMC approach 
than found in TQRPA model. Correspondingly the RPA 
calculation on top of the SMMC occupation numbers pre- 
dicts more GT strength in the energy range around 10 
MeV (corresponding to the OpOh centroid in the TQRPA 
calculation) which, as we will show below, results also 
in higher capture cross sections at low temperatures. In 
passing we note that, at low temperatures, the SMMC 
predicts a larger number of neutron excitations, but a 
smaller number of proton excitations across the pf-sdg 
shell gap, then the TQRPA model. This underlines the 
importance of pn-correlations, induced by isovector pair- 
ing and quadrupole interactions in the SMMC approach. 
Such pn-correlations are not considered in the present 
TQRPA calculations. 

In the SMMC the many-body correlations induced by 
the pairing+quadrupole interaction also yield a signifi- 
cantly smaller temperature dependence in the occupa- 
tion numbers than observed in the TQRPA approach. 
In this model the energies of the unblocked GT + tran- 
sitions essentially depend on temperature: at T > T cr , 
when unblocking is due to thermal excitations (thermal 
unblocking), they are smaller than the ones at T < T cr , 
when unblocking is due to configuration mixing. Obvi- 
ously, the significant shift of the GT + peaks to lower 
excitation energies favors electron capture. One can con- 
clude that at T > T cr GT + transitions in neutron-rich 
nuclei are more unblocked than at T < T CT . 

To explain the second effect, we consider the total 
strength for the j p — + j n transition 

St(jp -> jn) = $ *(JP ~> 3n) = (fjpjj^jpi 1 ~ n 3n)> 

(24) 



where ny T is the proton (r = p) or neutron (r = n) 
occupation factor 

nj = (0(T); qp|at m a. m |0(T); qp) = u)y) + v)x). (25) 

The value of nj (1 — nj n ) determines the unblocking 
probability for the j p — ► j n transition. As it follows 
from Eq. (|25[) the unblocking probability depends on the 
temperature through the coefficients of both the Bogoli- 
ubov and thermal transformation, i.e., it is determined 
by both configuration mixing and thermal excitations. 
We note again that at < T < T CT the total strength S t 
of the unblocked particle-particle or hole-hole transition 
is not concentrated in only one peak, but as it follows 
from (|2"3")) . is fragmented into four parts. Fig. EIe),(f) 
show the unblocking probabilities for the lfy 2 ~ * I/5/2 
and l<7g/ 2 ~~ * lff"/2 transitions as a function of tempera- 
ture. 

As seen from the figure the unblocking probabilities for 
both transitions have a minimum at the critical temper- 
ature. It is apparent that this minimum occurs because 
at T cr pairing correlations vanish while thermal effects 
are not yet sufficiently strong to occupy the lgg/2 pro- 
ton orbit or unblock the I/5/2 neutron orbit. As a result 
the total transition strength St decreases in the vicin- 
ity of the critical temperature. In contrast, this min- 
imum is absent in the SMMC unblocking probabilities 
(see Fig. |5[e),(f)). Here the residual interaction intro- 
duces a slight, but gradual increase of the probability 
with temperature. At T > 1.5 MeV the SMMC and 
TQRPA results converge as is expected in the high tem- 
perature limit. 

The fact that crossing shell gaps by correlations is a 
rather slowly converging process which requires the con- 
sideration of multi-particle-multi-hole configurations has 
already been observed in large-scale diagonalization shell 
model calculations, e.g. studying the calcium isotope 
shifts 45] or the Ml strength in argon isotopes (46j . 

For neutron-rich nuclei the contribution of first for- 
bidden p — > n transitions to electron capture is not 
negligible [3, [2(|. The strength distributions of first- 
forbidden 0~, 1~, and 2~ transitions in 76 Ge are shown 
in Fig.[7]for temperatures T = 0.2 and 1.3 MeV. The dis- 
tributions have been folded by the same procedure used 
above for the allowed transitions. As is seen from the 
figure, a temperature increase weakly affects the peaks 
in the _ , 2 _ strength distributions. The reason is that 
these are dominated by particle-hole transitions whose 
energy depends only weakly on temperature (in contrast 
to particle-particle and hole-hole transitions). With in- 
creasing temperature the peaks slightly shift to lower ex- 
citation energies due to the vanishing of the pairing cor- 
relations and some transition strength appears below the 
zero temperature threshold due to thermally unblocked 
transitions. 

Finite temperature induces a significant spread in the 
1~ transition strength distribution. The spread can be 
easily explained. At T — 0.2 MeV the main peak in the 
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FIG. 7: Folded strength distributions of first forbidden 0~, 1~ 
and 2~ p — > n transitions in 76 Ge at T = 0.2 MeV (left panels) 
and T = 1.3 MeV (right panels); _E is the transition energy. 
The strength distributions for the 2~ multipole correspond 
to 25 MeV electrons. St is the total strength. The arrows 
indicate the zero temperature threshold. The letters label 
the 1" transitions: A = l/£ /2 -> 2d™ /2 , B = 1/* -> 1$ 



C = I/7/2 



/7/21 



distribution is generated by three single-particle transi- 



tions: 



1 /7/2 



K/2> 1 / 5 /2 ^ l#/a. and 1/* 



7/2 



Iff 



9/2- 



The first is a particle-hole transition and its energy de- 
pends only slightly on temperature. The second and 
third ones are particle-particle and hole-hole transitions, 
respectively. As discussed above, the energies of particle- 
particle and hole-hole transitions at T > T cr are no- 
ticeably lower than that at T < T cr . Therefore, at 
T = 1.3 MeV, the peak is fragmented into three parts, 
resulting in a broadening of the 1~ strength distribution. 
The 1~ peak at E = 19 MeV is generated by the particle- 
hole transition 1/y / 2 ~* ^7/2 an( i) hence, its position and 
strength almost do not depend on the temperature. 

To reveal the importance of thermal unblocking for 
GT + transitions in neutron-rich nuclei we perform elec- 
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FIG. 8: (Color online) Electron capture cross sections (upper 
panels) for 76 > 80 > 80 Q e calculated within the TQRPA approach 
for various temperatures. Relative contributions of allowed 
transitions to the electron capture cross sections are shown in 
the lower panels. 



tron capture cross sections calculations. Within the 
present approach, the total cross section for capture of 
an electron with energy E e on a nucleus with charge Z 
is given by 



a(E e ,T) 



r 2 

2r7 



F(Z,E e )J2(E e 

Ji 



E 



(+h2$(+) 



Ji 



(26) 



where G w is the weak interaction coupling constant; 
F(Z,E e ) is the Fermi function that accounts for the 
Coulomb distortion of the electron wave function near 
the nucleus (see, for, example, @). Only allowed and 
first forbidden transitions are involved in the sum over J 
in the present study. 

In Fig. [5] the electron capture cross sections 

for 76,78,80 Ge 

are shown for three different tempera- 
tures. The temperature dependence of the cross sec- 
tions is most pronounced at moderate electron energies 
(Ee < 15 MeV): for E e — 15 MeV a temperature increase 
from 0.5 MeV to 1.3 MeV results in an enhancement of 
the cross sections by an order of magnitude. No such en- 
hancement was found in [2(| (see below). To make clear 
the reason of this enhancement, we calculate the rela- 
tive contribution of allowed transitions to the electron 
capture cross sections. The results are displayed in the 
lower part of Fig. [H 

It is seen that at E e < 15 MeV electron capture is 
mainly mediated by allowed transitions. Consequently 
the cross section enhancement is caused by thermal un- 
blocking of GT + transitions (the Fermi contribution to 
the cross sections is negligible). Furthermore, because 
of thermal unblocking, the electron energy below which 
electron capture is dominated by allowed transitions, 
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FIG. 9: (Color online) Electron capture rates for 76 > 78 ' 80 Ge calculated using the TQRPA approach as a function of temperature 
and for selected values of density pY e (in gem -3 ). For comparison, the rates obtained by the hybrid model are also shown. 



shifts to higher values: at T — 0.5 and 0.9 MeV this 
energy is 16 - 18 MeV, while at T — 1.3 MeV it is 
about 25 MeV. For larger electron energies, the first for- 
bidden transitions become increasingly important. As 
the strength of the first forbidden transitions is less sen- 
sitive to temperature change, the capture cross sections 
at E e ~ 30 MeV depend only weakly on temperature. 

The strong temperature sensitivity of the cross sec- 
tions at low electron energies reflects the temperature 
dependence of the TQRPA GT strength distribution, as 
discussed above. This has mainly two reasons. First, 
at low temperatures the dominant GT strength resides 
at larger excitation energies than in the hybrid model. 
Second, in the TQRPA the GT centroid shifts by several 
MeV, which is not observed in the hybrid model. Both 
facts, amplified by the strong phase space energy depen- 
dence, lead to a much stronger dependence of the cross 



section in the TQRPA model than in the hybrid model. 
As the GT contribution to the cross sections is larger 
in the hybrid model than in the present TQRPA calcu- 
lation, allowed transitions dominate to higher electron 
energies in the former. 

Fig. [9] compares the capture rates for 76 > 78 . 80 Ge as ob- 
tained in the hybrid and the TQRPA models. We note 
that the hybrid model rates are noticeably larger than 
the present ones at low temperatures. This is due to 
the increased unblocking probability in the hybrid model 
caused by many-body correlations which lead to a larger 
GT strength at lower excitation energies than in the 
TQRPA approach. With increasing temperature and 
density the differences between the rates of the two mod- 
els become smaller. This has two reasons. At first, with 
increasing temperature and density the average electron 
energy increases and the rates become less sensitive to 
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details of the GT strength distribution. Secondly, the 
GT strength distributions as calculated in the two mod- 
els become more similar with increasing temperature as 
discussed above. 

We note that the rates obtained in both models for the 
temperature and density regime, in which neutron-rich 
nuclei like those studied here dominate the composition 
during supernova core collapse (T > 1 MeV, p > 5 x 
10 10 g/cm 3 ), are large enough so that electron capture 
on nuclei dominates over capture on free protons as has 
been predicted in (20j . 

VI. CONCLUSION 

In this work we have considered the case of GT + and 
first-forbidden transitions in hot nuclei. We have applied 
the proton-neutron quasiparticlc RPA extended to finite 
temperature by the Thermo-Field-Dynamics formalism. 
The presented approach allows to treat charge- exchange 
transitions in nuclei at finite temperature without ap- 
plying Brink's hypothesis. It fulfills the Ikeda sum rule 
at finite temperature. As an example, we have calcu- 
lated the strength distribution for GT + transitions in 
54:56 Fe. We have observed the downward shift of the 
GT + strength with increasing temperature. This shift is 
caused by vanishing of pairing correlations and the ap- 
pearance of negative- and low-energy transitions. The 
shift of the GT + strength results in more enhanced elec- 
tron capture rates at high temperatures as compared to 
those obtained from shell-model calculations. We have 
found that the contribution of negative- and low-energy 
transition, i.e., transitions from thermally excited nuclear 
states, to electron capture is non-negligible even at low 
temperatures. We have also calculated the GT + strength 
distribution in the neutron-rich 76 ,78,80Q e nuc i e j_ was 
found that the temperature increase leads to a consider- 
able (of the order of 8 MeV) downward shift of the peaks 
in the strength distribution and reduces the total transi- 
tion strength in the vicinity of the critical temperature. 
This makes the unblocking effect for neutron-rich nuclei 
quite sensitive to increasing temperature in our model 
which is clearly observed in the electron capture cross 
sections and rates for 76,78 ' 80 Ge. 

If we compare our results to those obtained within 
the diagonalization shell model for the iron isotopes and 
within the hybrid SMMC+RPA model for the neutron- 
rich germanium isotopes, the importance of many-body 
correlations beyond those induced by pairing in our 
TQRPA model become apparent. For the isotopes 54 - 56 Fe 
Pauli unblocking is unimportant as GT transitions are 
possible even in the Independent Particle Model without 
correlations. The TQRPA describes the centroid of the 
GT strength rather well. However, it misses the low-lying 
GT strength, which is induced by multi-nucleon correla- 
tions. Such low-lying strength in 54,56 Fe, which is exper- 



imentally observed and reproduced by the shell model, 
is important for calculations of electron capture rates at 
low temperatures. Its neglect leads to underestimation 
of the rates. 

Pauli unblocking is crucial for the calculation of the 
GT strength and the associated electron captures rates 
for the neutron-rich germanium isotopes. Previous di- 
agonalization shell model calculations have shown that 
such cross-shell effects are rather slowly converging with 
increasing correlations across the shell gap and require 
the consideration of many-nucleon correlations. This is 
in line with the observation that the SMMC approach, 
which accounts for many-body configuration mixing, re- 
covers significantly more excitations across the pf—sdg 
shell gap than is found within the TQRPA approach, 
which, at low temperatures, derives Pauli unblocking 
mainly from 2p2h pairing correlations. These differences 
reflect themselves in different GT distributions and cap- 
ture rates at low temperatures and densities. However, 
the two models predict rather similar capture rates for 
the collapse conditions at which Pauli unblocking mat- 
ters, making the capture on nuclei dominate over capture 
on free protons. 

In summary, we have presented here a method which 
allows to calculate stellar electron capture rates at fi- 
nite temperature in a thermodynamically consistent way. 
This virtue makes it conceptually superior to the hy- 
brid approach of SMMC+RPA which has previously been 
used to estimate such rates for neutron-rich nuclei. In 
the present application of the model we have done a 
first step towards its complete realization describing cor- 
relations within the TQRPA. While this already recov- 
ers much of the essential physics, the detailed compari- 
son to shell model results implies that the model must 
be extended to include correlations beyond TQRPA. 
This can be achieved by taking into account coupling 
with complex thermal (e.g., two-phonon) configurations. 
For charge-exchange excitations in cold nuclei this prob- 
lem was resolved within the quasiparticle-phonon nuclear 
model [30l . l3l| and other approaches (J?} ■ It was found 
that the coupling with complex configuration strongly af- 
fects the RPA strength distribution. As discussed above, 
the details of the GT strength distribution are of particu- 
lar importance for weak interaction rates at low densities 
and temperatures. It is also desirable to extend the ap- 
proach to more microscopic effective interactions and also 
to consider the case of deformed nuclei. 
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